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We propose a new technique for determining the spin of new massive particles that might be 
discovered at the Large Hadron Collider. The method relies on pair-production of the new particles 
in a kinematic regime where the vector boson fusion production of color singlets is enhanced. For 
this regime, we show that the distribution of the leading jets as a function of their relative azimuthal 
angle can be used to distinguish spin-0 from spin-| particles. We illustrate this effect by considering 
the particular cases of (i) strongly-interacting, stable particles and (ii) supersymmetric particles 
carrying color charge. We argue that this method should be applicable in a wide range of new 
physics scenarios. 



I. INTRODUCTION 

The TeV energy scale is expected to contain an array 
of new particles predicted in extensions of the Standard 
Model (SM) that address electroweak symmetry break- 
ing, naturalness, and hierarchy problems. Among the 
most widely discussed scenarios are supersymmetry, ex- 
tra dimensions, grand unified theories, and technicolor. 
Assuming that new particles are discovered, the chal- 
lenge will be to ascertain which new physics model - 
including the possibility of one not yet considered - best 
accounts for them. This will be complicated by a de- 
generacy among many scenarios in basic experimental 
signatures. A determination of the particles' quantum 
numbers will provide key input in discriminating among 
the candidate models. 

One of the most important measurements - and the 
focus of this paper - is the determination of particle 
spin. It is, for example, one of the primary ways to dis- 
tinguish supersymmetry from universal extra dimensions 
(see e.g. PQ). Theorists have devoted considerable effort 
to developing methods for diagnosing spin at a hadron 
collider or a possible new linear collider 2 8 . The sim- 
plest method is to determine the kinematic dependence of 
the Drell-Yan pair production cross section [4] . More so- 
phisticated methods rely on the spin-dependence encoded 
in the decays of the produced particles. This approach 
is challenging for pp collisions at the Large Hadron Col- 
lider (LHC) since a full kinematic reconstruction of the 
decay is not possible in events with two or more invisi- 
ble particles. As many extensions of the SM contain dark 
matter candidates that would be produced in cascades at 
the LHC, this is an especially limiting constraint. Conse- 
quently, one must relate the spin-dependent correlation 
of daughter particles in the decaying particle rest frame 
to an appropriate quantity in the lab frame. The en- 
suing loss of kinematic information renders this method 
less than fully general, though it can be quite powerful in 
certain circumstances [3] [SHI] • Similar considerations 
apply to the measurement of spin information contained 
in quantum interference between different decay ampli- 
tudes [5]. It is therefore desirable to develop complemen- 



tary spin probes free of these limitations. 

In this paper, we propose a new spin measurement 
technique that does not require even partial reconstruc- 
tion of the new particles' momenta. The method relies 
on isolating a kinematic regime in which vector boson fu- 
sion (VBF) pair production is enhanced relative to other 
mechanisms (e.g., Drell-Yan). In this regime, we focus 
on the associated forward (or leading) jets and show that 
the differential distribution as a function relative jet-jet 
azimuthal angle, A0, can be used to diagnose the spin- 
statistics of the produced massive particles. In partic- 
ular, we find that the sign of the cos2A0 term in the 
distribution differs for spin-0 and spin-1/2 particles. The 
dependence on A0 can be understood as the physical 
observable arising from the interference of helicity am- 
plitudes associated with the vector bosons involved in 
the pair production of new particles. Information on the 
new particles' spin - encoded in the Lorentz structure 
of their pair production amplitude - is communicated to 
the interfering vector boson helicity amplitudes. The re- 
sulting azimuthal distribution of the leading jets probes 
this information. As we show below, by experimentally 
determining the sign of the cos 2A</> component of the 
dijet distribution, one has as a result particularly clean 
"diagnostic" of the new particle's spin. 

Our investigation of this technique has been heavily 
influenced by previous studies of dijet correlations in 
VBF [9-13 . The earliest work concentrated on the use 
of weak boson fusion (WBF) to discover an invisibly de- 
caying Higgs boson or to probe the CP properties of the 
Higgs coupling to the weak bosons [HI [TO] • At the par- 
tonic level, the WBF process can be isolated from the 
large QCD and tt dijet backgrounds using the distinctive 
pscudorapidity difference (A77) and dijet invariant mass 
(nijj) distributions. The resulting signal associated with 
the invisibly decaying Higgs involves an enhancement of 
the differential cross section for Ad> < tt/2, whereas the 
remaining Standard Model dijet background favors the 
opposite hemisphere. Subsequent Monte studies by the 
ATLAS collaboration indicate a persistence of this signal 
after accounting for parton showering and various trig- 
ger options [14]. A discovery of the invisibly decaying 
Higgs would be possible with 30 fb _1 of integrated lu- 
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minosity for a wide range or Higgs boson masses if the 
effective cross section is more than 60% of the SM Higgs 
production cross section. 

Subsequent studies of the Higgs boson CP properties 
considered both WBF [10] and gluon fusion (GF) [TUCE]. 
The Ari and nijj distributions for the two classes of VBF 
processes are quite distinct, leading to different event se- 
lection criteria for each. For both processes, however, 
the appearance of a distinctive Acj) distribution associ- 
ated with various CP-even and CP-odd couplings of the 
Higgs to vector bosons appears to be associated with 
leading jets that are widely separated in pseudorapidity 
(see, e.g., Fig. 11 of Ref. [E])- Thus, imposing an event 
selection cut on | appears advantageous for enhancing 
the A</> signal, even for GF for which the Ar] distribution 
itself is not useful for isolating the fusion process. 

Motivated by these observations, we consider the A(j> 
distributions for pair production of new particles through 
the VBF process in the large Arj regime. We find that 
as the Arj cut is relaxed, the distinctive signal associated 
with spin-0 or spin-1/2 particles vanishes. To gain further 
insight into this result, we show analytically that the sig- 
nal is associated with the two fusion bosons in the color 
(or weak isospin) singlet state. In the case of GF, color 
singlet dominance - or color "coherence" of the leading 
jets - has been intuitively associated with a rapidity gap, 
since a color singlet object is unlikely to radiate signif- 
icantly in the central region (for a discussion of these 
expectations, see, e.g. Ref. [TS] and references therein). 
Overlapping minimum bias events and minijet activity 
will populate the central region - considerations that are 
likely to apply to the earlier GF studies [TTJ [T^] as well 
as in the present case. Nevertheless, it is intriguing that 
the spin-dependent signal associated with the underlying 
hard event is associated with both a large separation in 
pseudorapidity and a color singlet configuration. 

As a practical matter, one must be able to distinguish 
the leading jets from jets produced in the decays of the 
produced particle pair, particularly in the case of GF 
where the latter is produced strongly. In this respect, 
we note that the forward jets relevant for our observ- 
able are essentially initial state radiation (ISR) associ- 
ated with the hard event. Typically, ISR is an unde- 
sirable but unavoidable fact of life at hadronic collid- 
ers; ISR jets being easily confused with jets originating 
in decays of new physics. Considerable work has gone 
into separating these two classes of hadronic activity (see 
e.g. Ref. [Tlfl H7] . in which ISR must be distinguished 
from the decays of gluinos in order to make mass mea- 
surements, and Ref. [18] where jet kinematics are used 
to statistically identify ISR jets). These studies suggest 
that the prospects are promising for isolating the jets of 
interest to the spin determination from the decay prod- 
ucts. We defer a detailed study of jet identification to 
future work, and concentrate here on the physics of the 
underlying signal. 

In order to minimize complications associated with SM 
backgrounds, we first demonstrate the spin-dependent 



A<f> correlation with a simulation for scalar and spinor 
stable strongly-interacting particles, i.e. i?-hadrons [El 
|2"0] . This is an interesting illustrative scenario in its own 
right, as the spin of stable particles cannot be measured 
using the standard techniques (see Refs. [5TJ [25] for ex- 
ceptions). We also show, through explicit simulations, 
that the same effect is present for strongly interacting 
supersymmetric particles and for Standard Model ti+jj 
production. In all these cases, the relevant vector boson 
is the gluon. In fact, ti + jj will be one of the primary 
backgrounds, as the corresponding forward jets will dis- 
play the same A<f> distribution as for jets associated with 
GF production of a new pair of massive fermions. It 
should be possible to minimize this background using a 
6-jet veto[E], a possibility that we will study in detail 
along with jet identification in a forthcoming publication. 
We anticipate that the same spin-dependent A(f> signal 
should also be present for weak boson fusion (WBF). Al- 
though WBF can be isolated through the Arj and m,jj 
distributions, and backgrounds further reduced through 
use of the central jet veto (see, e.g., Refs. [23H26] and ref- 
erences therein), the event rate would be very reduced. 
Consequently, we concentrate here on the GF case in or- 
der to demonstrate the basic concept of the technique. 

In the next Section we explain in detail the underly- 
ing concept behind our new spin measurement technique, 
and analytically demonstrate its utility in the toy ex- 
amples of scalars and spinors charged under an abelian 
gauge group. In Section |III[ we apply the technique to 
simples models with non-abelian QCD interactions, us- 
ing MadGraph/MadEvent [17] and CalcHEP simulations 
[2"5] . As will be shown, inclusion of color factors greatly 
complicates the analytic calculation, but does not change 
the overall applicability of the method. 



II. ANALYTIC CALCULATIONS WITH 
ABELIAN EXAMPLES 

As outlined in the Introduction, we are interested in 
the pair production via VBF of two heavy particles with 
two forward jets. The interference of the vector bosons 
helicities is set by the Lorentz structure of the pair pro- 
duction matrix element, and in turn affects the distribu- 
tion in A<fi of the two forward jets. To make this causal 
chain clear, we begin with an abelian gauge group as the 
vector bosons which fuse, and consider only the dominate 
Feynman diagrams in order to make the computation 
tractable. The full calculation of abelian gauge groups, 
kinematic cuts, and all involved diagrams requires simu- 



lation to perform, and will be described in Section III 



To set the theoretical framework, we consider pair pro- 
duction of massive particles in a regime where VBF is 
kinematically favored. As illustrated in Fig. [T] the event 
contains two forward jets (Ji and J 2) as well as the two 
new particles (each having mass m and momenta p\ and 
P2). The azimuthal angle <pj of each jet is defined rela- 
tive to a reference plane. The location of the reference 
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FIG. 1: A VBF event in an idealized detector volume. 
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FIG. 2: Kinematics for VBF in Frames I, II and X. Azimuthal 
angles <f>i and 02 are denned relative to the production plane 
of the massive particles (momenta pi and P2). Polar angles 
6\, 62 and 6 are defined relative to the beam axis in their 
respective frames. 



the contribution to A\ from the interfering amplitudes 
[9]. The fusion process tends to knock the produced 
pair and recoiling jets on opposite sides of a suitably de- 
fined reference plane, favoring a smaller A<p. In resonant 
production, this A<j) bias can be exploited in searches 
for "invisibly" decaying Higgs bosons produced through 
WBF [21 [T3]. In contrast, the cos2A0 term is dominated 
by the interference of the VBF amplitudes, receiving no 
significant kinematic contribution. Importantly, we find 
that the sign of the coefficient function A 3 is opposite for 
scalars (positive A 3 ) and spinors (negative A3). While in 
this Section, we only demonstrate the sign flip between 
the two spin choices, the results of our simulations in 
Section [TTT] will show that, while the cos2A^> correlation 
is subleading, it is sufficiently large to make its relative 
sign a useful spin diagnostic. 

The coefficient A3 arises from the interference of the 
VBF pair production matrix elements -M pa ir that differ 
by two unit of boson polarization: 

A 3 = (VS) Yl (Mi(h 1 )M 2 (h 2 )M piiil (h 1 ,h 2 )) 

\h i -h / . |=2 

( y M 1 (h , 1 )M 2 (h' 2 )M pail (h' 1 ,h , 2 )y . (2) 

Here (VS) is the phase space factor; ftAi and M. 2 are the 
parton — > parton 4- gauge boson production amplitudes 
(see Fig. |2|; and A^ pa ir is the VBF amplitude for pro- 
duction of the heavy pair. Performing the helicity sums 
leads to 



plane is arbitrary, as is direction of the z-axis, making 
the labeling of jets as "I" and "2" arbitrary as well (in 
the absence of jet tagging). Consequently, the only phys- 
ically quantity is the relative azimuthal angle A<p. 

Our labeling of the kinematics of the VBF amplitude 
is indicated in Fig. [2j largely following the conventions of 
Ref. [T3]. The polarization vectors of the gauge bosons 
are proportional to e *^*, where i = 1,2, and hi is the 
polarization of the boson defined relative to the posi- 
tive z-axis. The matrix elements for the emissions of 
the gauge bosons can be expressed as Mi(6i, 4>u hi) = 
e lh ''l> i Mi(8 l , 0, hi). Thus, any dependence on <fii and <p 2 
in the amplitude must be due to interference between 
matrix elements of the different polarizations. The re- 
sulting contribution to the differential cross section then 
yields set of 24 functions (sines and cosines) whose ar- 
guments are linear combinations of the <pj . As discussed 
above, in this paper we retain only the terms dependent 
on A(f>. Moreover, if the pair production amplitude is 
CP-invariant as we will assume here, then only the co- 
sine terms survive. Thus, 



da 
dA~4> 



A + Ai cos A(f> + A 3 cos 2A</> 



(1) 



The VBF kinematics induces a basic cos A0 depen- 
dence on the cross section that typically overshadows 



A3 = (VS)M l (+l)M 1 (-l)*M 2 (-l)M 2 (+l)* (3) 

X .Mp a ir(+1, -l)AV ir (-l, +1)* + (1 -!)• 

In the illustrative simple case considered in this Sec- 
tion, that of abelian gauge bosons (e.g., photons), the 
fusion amplitude A^ pa ir is generated by the standard 
Compton diagrams (with gauge coupling g). For spin- 
particles of charge Q, a straightforward computation 
gives 



M 



scalar 



g 2 Q 2 {(ti-e 2 )~ 

, (Pi ■ ^2) (pi - 92) ■ ei 



(Pi ■ e i)(Pi ~ Qi) ■ £2 



4 



2pi ■ qi 



4 



2pi ■ q 2 



(4) 



Choosing polarization vectors (with hi = — 1,0, +1) 
to reflect the angular momentum in the z-direction, 
ef = ef = e ± . It is then obvious that 
M.(+l,—l) = M(— 1, +1), which demonstrates that 
A 3 oc |A4 S caiar(il) =Fl) | 2 , and is therefore positive. 
Though a more lengthly calculation is involved, this ar- 
gument also holds for VBF scalar pair production by non- 
abelian gauge bosons. 

The calculation for spinors is similarly straightforward, 
but the origin of the negative sign of A 3 (opposite that 
of the scalars) is less transparent. However, a closed 
form expression can be obtained in the limit that the 
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two abelian gauge bosons are are ma ssless and o n-shell 
ill = Qi = 0). Defining A = \J\ — 4m 2 Js where 
s = (ki + k 2 ) 2 , we find 

. 64m 2 / 4m 2 , , \ . , 

A 3 oc 1 H — tanh A , (5) 

s \ sX J 

which is clearly negative. It is interesting to note 
that the minus sign (relative to the scalar calcu- 
lation), arises from the anticommutation of fermion 
operators in the simplification of the expression for 
■A^ pa i r (+1, — l)-M P air(— 1, +1)*, thus making the connec- 
tion between the sign of cos 2A</> and spin clear. 

We now briefly mention the difficulties in analytic cal- 
culation for the non-abelian fermionic example, before 
continuing on to the numeric simulations. In the ex- 
plicit calculation of the non-abelian example, two terms 
arise: a positive contribution with an antisymmetric color 
factor [T a ,T b ], and a negative contribution - identical 
to Eq. Q up to an overall a symmetric color factor 
{T a ,T b }. These two terms correspond to the two glu- 
ons being in a color octet and color singlet state, respec- 
tively. As we shall see, numeric calculation for off-shell 
gauge bosons confirms that the coefficient is overall neg- 
ative, meaning that the color singlet state (symmetric 
term) dominates. This is an intriguing result, and we 
give a heuristic explanation. Helicity conservation and 
the kinematics of nearly forward jets (typical of ISR) fa- 
vors a spatially symmetric configuration for the fusing 
gluons, and thus a color symmetric (singlet) two gluon 
state. However, without an application of kinematic cuts 
on the 4-body phase space, this is not reflected in the 
analytic calculation. 

III. SIMULATIONS WITH QCD GAUGE 
FORCES 

We therefore turn to a full calculation of the jet distri- 
butions using non-abelian gauge boson. As outlined pre- 
viously, a realistic computation involves additional sub- 
tleties, though we find from simulations that the sign 
difference persists. The non-abelian Compton amplitude 
A'fpair contains a dependence on the structure constants 
and generators of the gauge group as well as contributions 
from a diagram involving fusion of two gauge bosons to 
a third, followed by production of the heavy pair. In ad- 
dition, the coefficients A , A\ and A 3 receive additional 
contributions from heavy particle pair production via a 
bremsstrahlung gluon. If production was primarily via an 
on-shell intermediary (i.e., a resonance), the A(f> signal 
would be governed by the coupling of the intermediary to 
the two gauge bosons rather than by the spin of the pro- 
duced pair. However, such production is not dominate in 
the cases we consider here. 

To illustrate the spin-dependent sign of the cos 2A<fi, 
we first specialize to VBF production of i?-hadron pairs 
at the LHC (for an overview and references, see e.g., 
Ref. These stable, strongly interacting particles 



would have striking signatures at the LHC, either appear- 
ing as "slow muons" that are highly ionizing [29, 30], or as 
stopped tracks in the detector [3TJ [32]. The background is 
therefore limited to detector errors, which we will ignore 
for this study. i?-hadrons were originally proposed in su- 
persymmetric models with a gluino or squark as the light- 
est (or next-to-lightest) supersymmetric particle [33Tl39| . 
though other theories can provide similar particles (for 
an excellent overview, see Ref. [29 ). Null results from 
searches from ALEPH [40], CDF [41], and LEP2 [42], 
exclude particles of masses less than ~ 200 — 250 GeV, 
depending on the theoretical assumptions made, while 
the LHC should be able to find i?-hadrons up to a few 
TcV in mass [29l [43]. In addition, since standard spin 
measurements cannot be applied to stable particles, new 
techniques would be necessary in order to study the spin 
of i?-hadrons [22] . 

We assume some new Z2 quantum number (i?-parity in 
supersymmetry) that requires i?-hadrons to be pair pro- 
duced. As strongly interacting particles, this production 
can proceed via VBF mediated by gluons ( "gluon fusion," 
or "GF"). Matrix elements for two heavy fermions or 
scalars plus two jets from pp initial states are generated 
in MadGraph, and Monte Carlo simulations in MadE- 
vent |27j : such simulations correctly include all interfer- 
ence terms as well as production from a bremsstrahlung 
gluon. To isolate the kinematics in which the GF process 
is enhanced relative to Drell-Yan, we impose the follow- 
ing cuts: 

Vh ■ Vh < 0, M < 5, \r ln - r) h \ > 4.2 
PT, h > 30 GeV, p Td > 20 GeV, > 500 GeV (6) 

kfl-hadron| < 2.1, PT,i?.-hadron > 50 GeV. 

Here, r\j is the jet pseudo-rapiditiy; Mjj is the dijet in- 
variant mass; and px denotes the transverse momentum 
of a given particle. The jet cuts are necessary in or- 
der to reduce the non-VBF events [HI [25] , while the R- 
hadrons cuts are of lesser importance, and simply require 
the heavy states appear in the barrel of the detectors. 

We generate 40,000 events for the LHC at ^ = 10 TeV 
for simple, 500 GeV, scalar and Dirac spinor i?-hadrons, 
in which the heavy particles only couple to gluons. After 
all cuts the spinor cross sections are 33 fb (spinor) and 
21 fb (scalar). The differential cross sections da/dA<fi for 
the two models are shown in Fig. [3j along with the best 
fits to Eq. (l]). 

It has been demonstrated that the VBF-isolating jet 
cuts outlined in Eq. ^ tend to induce a positive cos A<f> 
mode in the distribution 9 . The expected positive A\ 
mode is found both the scalar and spinor simulations. 
The cos 2A<f) coefficient A 3 is not affected by the cuts, and 
we see a significant difference between the two spin cases. 
The scalar predicted, has a positive coefficient, 

while the spinor has a negative one. Appropriately nor- 
malized to the constant term, we find the scalar A 3 /Aq 
is 0.22, while in the spinor example, A 3 /A a — —0.14. 
Without background, the statistical significance can be 
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FIG. 3: Differential cross section da/dA(f> as a function of 
difference in the azimuthal angle of VBF jets A(j>. The differ- 
ential cross section for fermionic R-hadron pair production is 
shown in blue, and the best fit to Ao + Ai cos A<j)+A3 cos 2A0 
is shown with a solid line. The scalar cross section is in red, 
and the best fit is shown with a dotted line. 
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FIG. 4: Differential cross section da/dA<j> as a function of 
A(j>. The differential cross section for gluino pair production 
is shown in blue (scale on left), and the best fit with a solid 
line. The right-handed up-type squark cross section is in red 
(scale to right), and the best fit is shown with a dotted line. 



estimated as ~ |j43/v1o| V^N, necessitating ~ 30 fb -1 of 
luminosity for measurement at several sigma. We find 
that the dominant partonic subprocess involves gluon 
radiation from light quarks. We also observe that the 
A<j) signal disappears when the pseudo-rapidity cuts are 
removed, indicating the dependence of this method on 
kinematically selecting a region in which the color-singlet 
VBF process is relatively enhanced. 

The correlation between spin and the sign of A$ is 
independent of the charge conjugation of the produced 
fermions. In the case of supersymmetry, for example, 
one may wish to determine the spin of strongly interact- 
ing superpartners (squarks and gluinos) using GF. In this 
case the gluinos are Majorana fermions, in contrast to the 
Dirac R-hadron fermions explored above. The results of 
an analogous simulation for 550 GeV up-type squarks 
and 600 GeV gluinos are given in Fig. [3] Here we again 
see the sign change of the cos 2A(f> term when comparing 
squark (A 3 /A = 0.24) and gluino (A 3 /A = -0.09) pair 
production. In contrast to the .R-hadron case, isolating 
the signal associated with pair production of superpart- 
ners will require additional cuts to suppress Standard 
Model backgrounds. As our focus here is on the basic 
signal rather than on its isolation, we defer a detailed 
study of background reduction to a future analysis. 

To ensure that the effect is not an artifact of our choice 
of event generator, we have repeated the calculation for 
the dominant partonic sub-process using Calchep |28j 
and confirmed the MadGraph/MadEvent results. As a 
final check we have computed Standard Model tt+jj pro- 
duction in GF kinematics and find a negative sign for A3. 
Depending on the heavy particle mass m, tt + jj could 
thus become a significant background that would mimic 
the A<j) signal associated with new fermions. Implement- 
ing a 6-jet veto may be used to reduce this background 



IV. CONCLUSION 

In this work, we have demonstrated a proof of principle 
for the use of the sub-leading A<p distribution as a spin di- 
agnostic in a kinematic regime where color-singlet GF ap- 
pears to be kinematically enhanced and where the heavy 
particle pair production is dominated by their lowest or- 
der gauge interactions. We anticipate that the method 
will generalize to higher spin states and to pair produc- 
tion through weak vector boson fusion. Strong produc- 
tion of i?-hadrons provides the cleanest illustration as the 
signal has a large cross section and is essentially back- 
ground free. Other cases - such as VBF production of 
supersymmetric particles - in which the heavy particle 
decays produce additional jets and/or missing energy will 
require additional cuts for appropriate signal isolation. 
As indicated above, recent investigations in this direc- 
tion [Tni [T7] are encouraging for the prospects of appro- 
priate jet identification in the context of GF. A detailed 
study of these issues, background suppression, and other 
practical considerations will appear in forthcoming work. 
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